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STRESS-INTENSITY FACTORS FOR A THICK-WALLED 
CYLINDER CONTAINING AN ANNULAR IMBEDDED OR 
EXTERNAL OR INTERNAL SURFACE CRACK* 

by 

Ragip Erdol and F. Erdogan 
Lehigh IMiversity, Bethlehem, PA 


The paper considers the elastostatic axisymmetric problem 
for a long thick-walled cylinder containing a ring-shaped internal or 
edge crack. Using the standard transform technique the problem is 
formulated in terms of an integral equation which has a simple Cauchy 
kernel for the internal crack and a generalized Cauchy kernel for the 
edge crack as the dominant part. As examples the uniform axial load 
and the steady- state thermal stress problems have been solved and the 
related stress intensity factors have been calculated. Among other 
findings the results show that in the cylinder under uniform axial 
stress containing an internal crack the stress intensity factor at the 
inner tip is always greater than that at the outer tip for equal net 
ligament thicknesses and in the cylinder with an edge crack which is 
under a state of thermal stress the stress intensity factor is a de- 
creasing function of the crack depth, tending to zero as the crack 
depth approaches the wall thickness. 


1. INTRODUCTION 


In studying the fracture problem in pressure vessels, pipes, and 

other cylindrical containers , if the wall thickness relative to the 

mean radius is sufficiently small, the related crack problem may be 

solved by assuming the structure to be a shallow shell. In this case 

the problems which lend themselves to analytical treatment are those 

involving relatively short circumferential and axial through cracks 

in (infinitely) long cylindrical shells (for review and references see, 

for example, [1]). The problems of a through crack in a thick -walled 

cylinder and a part- through crack of finite size in a thick or thin- 
* 
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walled cylinder at the present time do not appear to be tractable. 
Nevertheless, these are the practical problems. Therefore, it would 
be of some interest if one can provide bounds for or give some idea 
about the magnitude and the trends of certain relevant quantities such 
as the stress intensity factors by solving the problem under certain 
idealized conditions regarding the crack geometry and the loading. One 
of the sin^lest of such idealized problems is considered in this paper. 
If (a) the flaw in the cylinder is a part-tlirough circumferential crack 
lying in a plane perpendicular to the axis, (b) the radial dimension of 
the flaw is relatively constant and its circumferential dimension is 
large compared to the wall thickness, and (c) the external loads are 
axisymmetric, then the axisyinmetric ring-shaped crack shown in Figure 1 
may approximate the problem around the center portion of the flaw. 

From the viewpoint of fracture propagation, since this mid-portion of 
the flaw is the most critical location, the approximate solution thus 
found would be very useful. The ring-shaped crack may be an internal 
crack (a<c<d<b) , an edge crack on the inner surface (a=c<d<b) , or an 
edge crack on the outer surface (a<c<d=b) (see Figure 1) . Once the 
governing integral equation is obtained, the problem can be solved for 
any axisymmetric quasistatic external load. In this paper the results 
will be given for uniform axial loading and for steady-state ten^er- 
ature distribution. 

The corresponding problem for the solid cylinder with a penny- 
shaped crack (i.e., a=0=c) was treated in [2-4], and for the strip 
with internal or edge cracks in [5] . 
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2. DERIVATION OF THE INTEGRAL EQUATIQN 

Consider the problem for a long thick-walled cylinder vdiicdi con- 
tains a concentric ring-shaped crack in z=0 plane, z being the axis of 
the cylinder (see Figure 1). Let the cylinder be subjected to (quasi- 
static) axisymmetric external loads which may be mechanical, thermal, 
or residual in nature and are independent of z. In the absence of the 
crack let the stress components at z=0 be 

a^ 2 (r, 0 ) = f(r) , o^^(r,0) = 0 , a<r<b , (1) 

v^ere f(r) is a known function. If the cylinder is sufficiently long 
the solution of the cylinder containing the crack may be obtained by 
superposition, namely, by adding to the (homogeneous) solution of the 
cylinder without the crack and under given external loads a perturba- 
tion solution obtained for the cracked cylinder in which self-equili- 
brating crack surface tractions are the only external loads. These 
tractions are equal and opposite to the stress given by (1) . From 
the viewpoint of fracture, the relevant problem is the latter. There- 
fore, only the perturbation problem will be considered in this paper. 

In terms of the Love function x(^>z) the axisymmetric problem 
may be formulated as follows [6] : 

V^x = 0 , (2) 

u(r, 2 ) = - , w(r,z) = [aCl-vJV^'x - §• ] 


C3) 



( 4 ) 



where u and w are r and z-components of the displacement vector, p is 
the shear modulus, and v is the Poisson's ratio. For an axisyinmetric 
elastic cylinder in which z=0 is a plane of symmetry and the stress 
state vanishes at z=+°° the biharmonic function y inay be expressed as 

X(r,z) =-[ [fx(s)I^(rs) +f 2 (s)rslj^(rs) 

■'o 

+ fjCs)K^(rs) + f^(s)rsK^(rs) ]sin zs ds 
f" -zp 

+ f,(p)p(2v+zp)e J (rp)dp , (5) 

■•o 

where the functions f^,...,f^ are unknown. From (3), (4), and (5) the 
displacements and the relevant stress con^onents are found to be 

-00 

u(r,z) = - ^ [f^Cs)Ij^(rs) + f 2 Cs)rsIj^(rs) 

- f 2 (s)K^(rs) - f^(s)rsK^(rs)]s^cos zs ds 

pOO 

f^(p)p^(l-2v-zp)e ^^J^(rp)dp) , (6) 

-00 

w(r,j) = ^ {| |^[f^(s)IgCrs) + f2(sH4Cl-v)I„(rs) 

+ rsl^(rs)] + f 2 Cs)K^(rs) + f^Cs) [rsK^(rs) 
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a (r,z) 


4(l-v)K^(rs) ]]s^sin zs ds 


-zp. 


f5(p)P [2(l-\))+zp]e J^(rp)dp} 


C7) 


2 

TT 


{s f^(s)[-sl.^(rs) + I^(rs)/r] 


+ s^f2Cs) [(2v-l)I^(rs) - rsl^(rs)] 


s^f^Cs) [sK^(rs) + K^(rs)/r] + s^f^(s) [(l-2v)K^(rs) 


•QO 

_ 3 . 


rsK^(rs)]}cos zs ds + J p f^Cp) [p(l-zp) J^Crp) 


-zp. 


(l-2v-zp)Jj(rp)/r]e dp 


( 8 ) 


2 

TT 


r 


{fj^Cs)I^(rs) +£2(s)[2(2-v)I^Crs) + rsl^(rs)] 


+ f2(s)K^(rs) + f^(s) [rsK^(rs) - 2(2-v)K^C^s) ]}cos zs ds 


r“ 


o 


5 Cp) Cl+zp) e (rp) dp , 


(9) 


2 

TT 


{£^(s)I^(rs) + £2(5) [rsl^(rs) + 2(l-v)I^Crs)] 


3 


- £2Cs)K^(rs) + £4(5) [-rsK^(rs) + 2(l-v)Kj^Crs}]}s sin zs ds 


c -zp 

p £5(p)ze J]^(rp)dp 


( 10 ) 
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In the crack problem shown in Figure 1 the following are the bound- 
ary conditions : 

CT^j.(a,z) =0 , a^^(a,z) = 0 , 0<z<“ , (11) 

cfj.j.(b,z) = 0 , ° 

cJrzC^'jO) = 0 , a<r<b , (13) 

= -fCr) , c<r<d , 

w(r,0) - 0 , a<r<c , d<r<b 

From (10) it may be observed that the condition (13) is satisfied by the 
solution (5) . The four homogeneous conditions (11) and (12) may be used 
to eliminate four of the unknowns f^ , . . . , f ^ and the mixed boundary con- 
ditions (14) may be used to obtain a system of dual integral equations 
for the fifth unknown function. Thus, using (8) and (10) the boundary 
conditions (11) may be expressed as 

f^(s) [-sl^(as) + I^(as)/a] + sf 2 (s) [(2v-l)I^(as) - asl^(as)] 

- f^(s) [sK^(as) + K^(as)/a] + sf^(s) [(l-2v)K^(as) - asK^(as)J 
1 r” P'^fsCp) 

= ^ “ (l-2v)J^(ap)/a 

t [pj^(ap) - J 3 ^(ap)/a] dp , (15) 



- 6 - 



f^(s)I^Cas) + f 2 (s) [asl^(as) + 2(l-v)I^(as) 


f 3 (s)Ki(as) + f^(s) [2(l-v)Kj^(as) - asK^(as}] 




2sp6fg(p) 

„ Ji(ap)dp 


( 16 ) 


Conditions (12) give two inore equations which may be obtained from (15) 
and (16) by replacing a by b. Now, rather than substituting 
given by these four equations into (14) and obtaining a system of dml 
integral equations for f^ (p) , the problem may be reduced to a singular 
integral equation in terms of a new unknown function defined by 


^w(r,0) = g(r) 


(17) 


From (7) , (17) , and the second equation o£ (14) it follows that 


P‘'f5Cp3 = I tg(t)J^(pt)dt . 


( 18 ) 


If we now solve (15) , (16) and two similar equations for r=b for 
f^,.^,f 4 in terms of g(t) by using (18) and substitute the results into 
the first equation of (14) , by using again (18) and the expression for 
*^7 7 given by (9), after somewhat lengthy but straightforward manipula* 
tions, we obtain 


rd 


+ k(r,t)]g(t)dt = - 7tf(r) 


c<r<d 


(19) 
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vdiere the kernel k(r,t) is given by 


k(r,t) = k^(r,t) + Ztk^Cr.t) 


k^C 


r t'l - m(r,t) 

t-r t+r 


fE(r/t) , r<t , 


m(r,t) =<, 


I 

jE(t/r) + 


tr 


K(t/r) , r>t 


ds ** 


^ B.h.)[2(2-v)I^(rs) 


+ rsl^(rs)] + (i Cj^h^)KQ(rs) + D^h^) [rsK^(rs) 


C20) 

( 21 ) 


( 22 ) 


- 2(2-v)K^(rs)]> . (23) 

Here, K and E are the conplete elliptic integrals of first and second 
kind, respectively and the functions A(s), A^(s), B^(s), C^(s), Dj^(s), 
and hj^(t,s), (i,l,..,4) are defined in the Appendix. The expressions 
of the infinite integrals used in this part of the analysis are given 
in [7], From the second equation of (14) and the definition (17) it 
is clear that the integral equation must be solved under the following 
single-valuedness condition: 

rd 

g(r)dr = 0 . (24) 

■'c 
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3. INTERNAL CRACK 


A close examination of the kernel k(r,t) defined by (20j to (23) 
would show that the first part k^^(r,t) has a simple logarithmic singu- 
larity (in the form of log|t-rl), whereas the second part k2(r,t) is 
bounded in the closed interval c<(r,t)<d provided a<c<d<b, that is if 
the crack is a fully imbedded internal crack. In this case the Cauchy 
kernel l/(t-r) is the dominant kernel, the index of the integral equa- 
tion is +1, the solution is of the form [8] 

g(r) = G(r) [(r-c)Cd-r)] ^ (25) 


and, hence a standard numerical technique such as that described in [9] 
may be used to determine the unknown function G(r) which is bounded in 
the closed interval c<(r,t)<d. Thus, defining the following normalized 
quantities 


d-c _ ^ d+c 


r = p + 


d-c . d+c 


, t = Y^T + 


g(r) = (j)(p) = 


F(pKl-p=^r’‘ , f(r) 


= P(p) , K(p,t) = 


(26) 


equations (19) and (24) may be expressed as 




1-v 


P(pj , 


p ""I 


j=l,..,n-l 


(27) 
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( 28 ) 


n 


I jF(Pi) =0 , 

i=l “ ^ 


- cos , X = 1,.. ,n 


Pj = cos ^ , j = 


(29) 


After solving the integral equation the stress intensity factors 
which are defined by 


k(c) = lim /2 (c-r) a (r,0) , 

r->c 


k(d) = lim /2(r-d) a (r,0) 


T-*d 


zz 


(30) 


may be obtained as follows: 


k(c) = lim /2(r-c) g(r) = ^ >/(d-c)/2 F(-l) , 1 


r-K: 


k(d) = -lim ^ >^Id^ g(r) = - ^ Ad-c)/2 F(l) • J 

r-»^ 


(31) 


whefe~’F{<iQ and F(+l) are obtained from F(xp, i=l,..,n, by using the 
interpolation fo^mtiias^iven in [10] . 


EDGE CRACKS 


In the case of edge cracks (i.e. » for a=c<d<b or a<c<d=b, Figure 1) 
the integral equation (19) is still valid. However, in this special 
case the kernel k(r,t) in (19) is no longer bounded in the corresponding 
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closed interval and, of course, the single-valuedness condition (24) is 
no longer valid. Therefore, the solution of the problem requires some 
care. The asymptotic analysis of the integrand in (23) for large values 
of s indicate that the kernel k 2 (r ,t) may be e:)q)ressed as the sum of 
two parts as follows: 

k2(r,t) = k2^(r,t) + k2g(r,t) (32) 


where k 2 £ is bounded in the corresponding closed interval and k 2 ^ becomes 
unbounded as r and t approach the end point c=a or d=b. After some manip- 
ulations the asymptotic expressions for the integrand and the singular 
part of the kernel k 2 are found to be 




o 


ds 

2>/ft 


[2s^(r-a) (t-a) - 3s(t-a) 


-s(r+t-2a) 

- s(r-a) + 2]e 


_ 1 r- 1 ^ 6(r-a) _ 4(r-a)^ ■■ 

*• r+t-2a (r+t-2a)^ (r+t-^a)*^' -• 


(33) 


in the neighborhood of the end point a (i.e., for r-a«d-a and t-a«d-a) , 
and 


kjsCr.t) = - 


[2s^(b-r) (b-t) - 3s(b-t) - s(b-r) 


0 2/rt 

-s(2b-r-t) 

+ 2]e ds 


= 1 r 1 _ 6(b-r) ^ 4(b-r)^ , 

2/rt 2b-r-t (2b-r-t)^ (2b-r-t) 


(34) 
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in the neighborhood of the end point b (i.e. , for b-r«b-c and b-t«b-G) . 
It should be noted that the singular kernels (33) and (34) arc essentially 
the same as the expressions found for the corresponding plane strain prob- 
lem considered in [5]. Thus, repeating the analysis of [5], it may be 
shown that the solution of the integral equation (19) , which now has a 
generalized Cauchy kernel, has no power or logarithmic singularity at the 
end point which is on the surface and the solution is of the following 
form 


g(r) = q(r)(d-r) 
g(r) = C^r)(r-c) ^ 


a=c<r<d<b , 
a<c<r<d=b , 


(35) 


where Giand G 2 are bounded in their respective closed intervals. In this 
case the stress intensity factors are given by 

k(d) = -im /zWT) g(r) = - ^ GjCd) , a=c , 

C36) 

k(c) = lira /2(c-r) gCr) = ^ G^Cc) , d=b . 
r->-c 

The integral equations for the edge cracks may be solved by follow- 
ing the procedure outlined in [5]. 

5. RESULTS AND DISCUSSION 

As an exanple, the integral equation (19) is solved under two types 
of input functions, namely, a uniform stress in z direction and a steady 
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state thermal stress distribution. In the former case 

= £(r) = 0 ^ , ,a<r<b (37) 

and in the latter (see, for example, [11]) 

0^[l-2 log (b/r) " log(b/a)]/log(b/a) (38) 

where 

0 ^ = aE(T^-T^)/2(l-v) , (39) 

a, E, T^, and being the coefficient of thermal expansion, Young’s 
modulus, outer wall temperature, and inner wall tenqierature , respect- 
ively. The calculated results are shown in Tables 1-4. Figures 2 and 
3 show some sample results which are taken from Tables 3 and 4. 

Table 1 shows the normalized stress intensity factor 

k' = k/(0^/(d-c)/2) (40) 

for a thick-walled cylinder which contains a ring-shaped internal crack 
along c<r<d and is under uniform axial tension ^ 2 z~^o 
may be seen that as the crack size becomes very small in comparison with 
other dimensions k’ approaches 1 which is the value for a uniformly 
loaded infinite plane containing a line crack. This was verified by 
substituting (d-c)/b=10"^ in the program. It may also be observed that 
if the thicknesses of the net ligaments, c-a and b-d are equal then the 
stress intensity factor k(c) at the inner tip is always greater than 
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k(d) at the outer tip, whereas in the corresponding plane problem for a 
strip they are, of course, equal. The practical consequence of this is 
that generally one would expect the internal flaws in thick -walled 
cylinders to propagate inward. Another observation which may be worth 
mentioning is that, for a fixed c/b, as d/b increases (starting with 
d=c), initially there is a slig^it reduction in k’(d) , then it increases 
and becomes unbounded as d->b. However, as seen from (40), the stress 
intensity factor k(d) itself is proportional to (d-c) ^ and, despite the 
slight decrease in k'(d), is always an increasing function of the crack 
size. If b=«> with a and c finite one would expect that as d increases 
k(d) would approach — which is the stress intensity factor for a 
penny- shaped crack in a uniformly loaded infinite medium. This means 
that in this case 2 /tt is the asymptotic value of k’(d) and l>k’ (d)>2/7T. 
For d-c«b the reduction in k’(d) is therefore not surprising. 

Table 2 shows the normalized stress intensity factors in a thick- 

walled cylinder with an internal crack which is under a state of steady- 

state thermal stress given by (38) and (39) . Here it is assumed that 

T^>T^, i.e., the cylinder is heated from outside. For again a fixed c, 

it is seen that as d increases k* (d) decreases and eventually becomes 

zero as d advances into the compressive stress region. Also as d-K:, 

k(d) and k(c) approach the local plane strain value i.e., k(d)->f(c) 

/(d-c)/2, k(c)^f(c)/(d-c)/2 where f(c) is the magnitude of the thermal 

stress a at r=c (see equations 38 and 39) . This was again verified 
z z 

by selecting (d-c)/b=10 
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Table 1. Stress intensity factors in a thick-walled cylinder with an 
internal crack under uniform axial stress 
></(,d-cJ/2 ) , k'Cc) = k(c)/Caov^(d-c)/2} . 


0^2=00, k'(d)=kCd)/(g^ 


d/b 

a/b=0.3, c/b=0.4 

a/b=0.3, 

c/b=0 . 5 

a/b=0.2. 

c/b=0 . 4 

k’(d) 

k ' (c) 

k'(d) 

k’ (c) 

k’(d) 

k’(c) 

0.4 

0.5 

-^1.0 

0.9973 

-s-1.0 

1.0591 

-^1.0 

-vl.O 

-^1.0 

0.9864 

-^1.0 

1.0405 

0.6 

1.0167 

1.1493 

0.9931 

1.0377 

0.9963 

1.1017 

0.7 

1.0626 

1.2690 

1.0150 

1.1035 

1.0345 

1.1900 

0.8 

1.1645 

1.4416 

1.0833 

1.2098 

1.1262 

1.3252 

0.9 

1.4361 

1.7374 

1.2870 

1.4022 

1.3726 

1.5644 

1.0 

-+O0 


->00 


-KX) 



Table 2. Stress intensity factors in a thick-walled cylinder with an 

intemal crack under steady state thermal stresses. k^(d)/(a^ 
y^Wc)[l ) , k{(c)/(at/(d-c)/2 ) , at=aE(To-Ti)/2(l-v), a/b=0.37 
c/b=0.4. 


a/b 

0.4 

0.401 

0.5 

0.6 

0.7 

0.8 


-^0.8893 

0.8859 

0.6044 

0.3762 

0.1827 

0.0093 


-J-0.8893 

0.8885 

...■ 

0.8328 

0.7836 

0.7279 

0.6591 


Table 3 shows the normalized stress intensity factors for an edge 

crack on the inside wall of the cylinder (i.e., a=c, d<b) . Here k'(d) 

refers to the uniform axial loading, k^(d) refers to the case of thermal 

stress, k' (d) is also shown in Figure 2 in order to conpare it with 

plane strain and penny-shaped crack results. It is seen that as d-m 

k(d)->1.1215 a 7d-a for uniform axial tension a =a and k(d)->1.1215 
o zz o 

f(a)/d-a for thermal stress f(r) given by (38). Note that for a half 


The limiting values of the stress intensity factors for the edge cracks 
given in Tables 3 and 4 were also verified by letting the crack depth 
to be 10 “^b, where a/b is of order unity. 
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Table 3. 


Stress intensity factors for the inside edge crack in a thick- 
walled cylinder (a=c<d<b, Figure 1). k* (d)=k(d)/ao /d-a (ex- 
ternal load: uniform tension cr 2 z=cJo) , k{ (d) =k (d) /ot /d- a (ex- 
ternal load: steady-state thermal stress, crt"^E(To~Ti)/2(l-v)) . 


— 

d/b 

a/b= 

3.3 

a/b=0.8 

lc’(d) 

k-(d) 

k^Cd) 

0.30 
0.40 
0.50 
0.60 
0.70 
0.80 
0.84 
0.90 
0.95 
1.00 
1 

->-1.1215 

0.9655 

0.9162 

0.9204 

0.9764 

1.1172 

1.4981 

-KO 

^1.5333 

1.0120 

0.7099 

0.5027 

0.3498 

0.2323 

0.1372 

->■0 .0 

->-1.2046 

0.9671 

0.7190 

0.4751 

->-0.0 


Table 4. Stress intensity factors for the outside edge crack in a thick- 
walle d cylinder, a/b=0.8 (a<c<d=b, Figure 1). k’ (c)=k(c)/ao 
yb^c (for uniform axial stress , k^ (c) =k (c) (for 

thermal stress, o-t-=aE(Ti-TQ)/2(l-v)) . 


c/b 

1.00 

0.98 

0.96 

0.94 

0.92 

mi 

0.85 

0.80 

1^ 

->-1.1215 

1.1615 

1.2513 

1.3706 

1.5141. 

1.6817 


-XX3 

n 

->-1.0383 

0.9487 

0.8905 

0.8385 

0.7832 

0.7189 

0.4892 

->-0.0 


plane with an edge crack of depth Z under uniform tension the stress 
intensity factor is 1.1215 c^/£ (see, for example, [5]). Thus the ex- 
pected results are recovered in limit. Table 3 and Figure 2 show that 
for a fixed a/b as d/b increases starting with a/b, the initial reduction 
in k* (d) is much more noticeable than the internal crack case. Figure 2 
also shows the stress intensity ratio for the penny-shaped crack in an 
infinite solid, namely 2/ir. The dashed line shown in the figure approx- 
imates the stress intensity factor ratio k' (d)=k(d)/a^»4i-a in an infinite 
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solid containing a cylindrical hole and an internal edge crack. It would 
be expected that for d-a«a^k’ (d)->-1.1215f for d»a^k' (d)->2/'n‘/ and, in be- 
tween ^k ’(d) would be a monotonically decreasing function of d. This 
appears to be the reason for the decrease in k’ (d) when d-a is siiiall 
compared to b. 

The result for the outer edge crack in a uniformly loaded cylinder 
is shown in Table 4. Here the stress intensity ratio k'(c) increases 
quite rapidly with increasing crack depth b-c. 

For a/b=0.8 the edge crack results for the steady -state thermal 
stress case are shown in Tables 3 and 4. The same results and the ther- 
mal stress a /a^ given by (38) are also shown in Figure 3. It is assumed 
that T >T. for the inner edge crack and T <T. for the outer edge crack, 

T^ and T^ being the outer and inner wall temperatures (Figure 3 shows 
only Tq>T^) . Note that (a) the stress intensity factor is always posi- 
tive if the edge crack starts from the wall which is under tension, (b) 
it monotonically decreases as the crack depth increases, and (c) it 
approaches zero as the thickness of the net ligament goes to zero. 
Basically, this is the same result found in [12] for a strip under (self- 
equilibrating) iesidual stresses. 
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APPENDIX 


The functions A^(s) , B^(s), Cj|^Cs), Dj,(s), h^(t,s), (i=l,.., 4 ), 
and A (s) . 

A^ = (-M^L2 + ’ 

*2 = • 

= Mj^[-sLj^Z^ •♦• L2 (sZ 2 + Z^^/a)] 

+ M2[-L4(sZ2 + Zj/a) + sZ^L^] - AK^(bs)/Z^ , 

\ = M^[L^(sZ3-Z^/b) - SZ4L2] 

+ M2[sL 4Z^ - L3CSZ3 - Z^/b) ] + AK^(as)/Z^ ; 

“ '^2^1 * 

B3 = “sL^Zg + ^2(322 + Z^/a) , 

B^ “ L2(sZ 3 - Zj^/b) “ sL2Z^ ", 

Cj^ = (-L2M3 + L^M^)Z^ , 

C2 = (L^M3-L3M4)Z^ , 

C3 = M3[-sLj^Z3 + 1^2 ^^^2 * ^ 

+ M4[-L^(sZ 2+ Z^/a) +SL3Z3] - AIj(bs)/Z^ , 
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= M3[L^(sZj-Zj^/b) - sL^Z^] 

+ M^[sL^Z^ - LjCsZj - Zj^/b) ] + fil3(as)/Z3 ; 

”l ' Vl ’ 

®2 ° '^3^1 ’ 

D3 — “L^(sZ2 ^ 3 3/1) + 5L3Z3 f 

” ^^4^4 ~ ^3^^^3 ” j 

Z3 = K^CasJI^tbs) - K3(bs)l3Cas) , 

Z2 = I^jCasjKjCbs) * KgCas)l3(bs) , 

Z3 = I^(bs)K3(as) + K^CbslI^Cas) , 

Z^ = I^CasJK^Cas) + K^(as)l3(as) , 

Zj = Ig(bs)K3Cbs) + Kjj(bs)l3(bs) ; 

Mj^ = [asI^(as)K^(bs) - sblQ(bs)Kj^Cas) - 2 (l-v)Z^]Z-^ , 

M2 = [-asK^(as)K^Cbs) + bsK^Cbs)K^(as) , 

= [asI^Cas)I^(bs) - bsI^(bs)I^(as)]Z'^ , 

^ [2(l-v)Z^ - asK^(as)I^(bs) +bsK^Cbs)I^Cas)]Z"^ j 
= as2[Z2K^(as) -Z^K^(as)] - bs^Z^K^Cbs) - Z^K^(as) » 
L2 = as^Z^K^Cas) - bs2[Z3K^(bs) + Zj^K^(bs)] - Z^K^Cbs) , 
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Lj = as^l-ZjI^Cas) -Z3^Ij^(as)]+bs^Z^I^Cbs) - Z, 

- -as^ZjI^(as) +bs^[Z3lg(bs) - Zj^I^Cbs)] - Z 

^ "" ^ 1^4 ■ ^ 2^3 ’ 

h^ = -5[I^(as)K^(ts)(2^^-as^) + s\l^(as)K^(ts) 

+ sI^(as)K^(ts) - f- Ij^(as)K^,Cts)] , 

hj = -s[I^(ts)lC^Cbs)(^i^-s2b) + ts^IgCts)KgCbs) 

- sI^(ts)KgCbs) + f- Ij,(ts)Kj(bs)] , 
hj = -s^[aI^(as)Kj^(ts) - tIj^[as)K^(ts) ] , 

= -s^[-bIj^(ts)K^(bs) * tl^Cts)Kj^(bs)] . 


I^Cas) , 
I^Cbs) ; 
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Figure 2 . The stress intensity factor for an (internal) edge crack in a 
thick-walled cylinder mder uniform axial stress 022=^0. 
a/b= 0 . 3 , k' (d)=k(d)/ag/d-a (the dashed line is the approxi- 
mate k’ (d) for an infinite medium with a cylindrical hole of 
radius a and an internal edge crack of depth d-a) . 





Figure 3. Stress intensity factors for an internal edge crack, k^(d) , 
and for an external edge crack, k^(c) , in a thick-walled 
cylinde r und er steady- stat e the nnal stress. a/b=0.8, k^(d)= 
k(d)/o^/d-a, k* (c)=kCc)/o^/b-c, o-(^=aElTo-TjLl /2(l-v] . 



